Recovering the superconducting state of superconductor-ferromagnet multilayer 
structures by spin accumulation above the Pauli pair-breaking magnetic field 
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We study theoretically the simultaneous influence of spin accumulation potential eVf — eVj. and 
the Zeeman exchange field on singlet superconductivity. It is shown that the pair-breaking effect 
of the Zeeman field can be fully compensated by creation of the appropriate spin accumulation 
potential in the superconductor. Moreover, superconductivity can be recovered for exchange fields 
well exceeding the Pauli limiting field. It is proposed that the effect can be experimentally realized 
on the basis of voltage biased junction consisting of a thin superconducting film sandwiched between 
two half metals. 
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One of the mechanisms destroying singlet supercon- 
ductivity is the Zeeman interaction of electron spins with 
magnetic field. The behavior of a magnetic superconduc- 
tor with an exchange field h was studied long agoi~— . It 
was found that homogeneous superconducting state be- 
comes energetically unfavorable above the paramagnetic 
(Pauli) limit h = A/y/2. As it was predicted 1 ^, in a 
narrow region of exchange fields exceeding this value su- 
perconductivity can appear as an inhomogeneous state 
with a spatially modulated Cooper pair wave function 
(LOFF-state). Thus, the exchange field is destructive to 
singlet superconductivity. An exception, proposed in the 
literature, is a special type of clean superconducting mul- 
tilayered system, where the paramagnetic limit can be 
enhanced and a nonuniform superconducting state can 
be induced under in-plane magnetic field 5 . 

Superconductor/ferromagnet (S/F) hybrid structures 
also can behave analogous to magnetic superconductors. 
In particular, it was shown^ that a thin S/F bilayer with 
thicknesses obeying the conditions dp -C £p and ds <C £s 
is equivalent to a magnetic superconductor with an effec- 
tive exchange field h e ff = hNpdp / '{Npdp + N$ds) and 
an effective superconducting order parameter A e ff 
AN s d s /(N F d F + N s d s ). 



Here £f = \/Dp/h and 



y/Ds/\A\ are the magnetic and superconducting 
coherence lengths, |A| is the superconducting order pa- 
rameter in the bulk material, Np,s denote the densities 
of states at the Fermi level for the ferromagnet and su- 
perconductor and Dp,s are the corresponding diffusion 
constants. Another way to "apply" an exchange field on 
a thin superconducting film is to contact it to a ferromag- 
netic insulator— ~—. It was observed experimentally^ and 
justified theoretically- 1 ^ that the effective exchange field 
induced in the film scales with dg 1 . 

The simultaneous applying of the exchange field and 
creation of spin-dependent quasiparticle distribution in 
S/F heterostructures can lead to qualitatively new phe- 
nomena. In particular, it was shown recentl y 12 ' 13 that 
creation of spin-dependent quasiparticle distribution in 
the interlayer of S/F/S junction (or S/N/S junction with 
magnetic S/N interfaces) leads to appearence of an addi- 
tional contribution to the Josephson current through the 



junction, which under certain conditions can enhance it 
considerably. In the present paper we demonstrate that 
for a thin superconducting film the destructive effect of 
the exchange field can be fully compensated by the cre- 
ation of spin-dependent quasiparticle distribution in it. 
This effect takes place even if the exchange field exceeds 
the paramagnetic limit considerably, that is under the 
condition that superconductivity of the equilibrium film 
is fully suppressed. 

As an example of a system, where the uniform ex- 
change field and the spin-dependent quasiparticle distri- 
bution can be realized simultaneously, we propose here 
a voltage-biased half metal/superconductor/half metal 
(HM/S/HM) heterostructure. A thin film (d s < £ s ) 
made of a dirty s-wave superconducting material is sand- 
wiched between two half-metallic layers with opposite 
directions of magnetization. Half-metallic behavior has 
been reported in CrO ^ 14 ' 15 and in certain manganites^i. 
In- plane effective uniform exchange field h e f f in the film 
is supposed to be created by spin-active interfaces with 
half metals. The spin-dependent quasiparticle distribu- 
tion in the film can be generated by applying a voltage 
bias between the two half metals. In this case for spin- 
up subband the main voltage drop occurs at one of the 
HM/S interfaces, while for spin-down subband - at the 
other. As a result, the distribution functions for spin-up 
and spin-down electrons in the superconducting film are 
to be close to the equilibrium form with different electro- 
chemical potentials. 

As we consider a non-equilibrium system, we make 
use of Keldysh framework of the quasiclassical theory, 
where the fundamental quantity is the momentum av- 
erage of the quasiclassical Green's function g(x, e) — 
(g(pf, x, e)) Pf . Here x is the coordinate normal to S/HM 
interfaces and e is a quasiparticle energy. The Green's 
function is a 8 x 8 matrix form in the product space of 
Keldysh, particle-hole and spin variables. In the super- 
conducting film it obeys the Usadel equation 

—d x (gd x g) + [eT 3 a p - AT 1 ia 2 pQ,g] = , (1) 
where r^, cr, and pi are Pauli matrices in particle- hole, 
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spin and Keldysh spaces, respectively. To, cto and Po 
stand for the corresponding identity matrices. Eq. ([1]) 
should be supplied with the normalization condition g 2 = 
— 7t 2 to<to/9o- -^ 1S convenient to express Keldysh part of 
the full Green's function via the retarded and advanced 
components and the distribution function: g K — g R (p — 
ipg A . The distribution function is diagonal in particle- 
hole space: if = tp(To + T 3 )/2 + <j2<P&2{To — Ts)/2. The hole 
component if of the distribution function is connected to 
if by general symmetry relation^ <f(s) = —<Ji(p{—E)o?,. 

We consider the case when there is the only magneti- 
zation axis in the system (the magnetization directions 
of the PM's are antiparallel) . Then there are no equal- 
spin triplet superconducting correlations in the system 
and all the matrices in spin space can be represented as 
sums of two spin subband contributions (a =t, ■!)• F° r 
later use we define here the anomalous Green's function 
fR ,A = (ff> A ( ao + CT3 )/ 2 + /«>Vo - ^)/2)ia 2 and the 
distribution function ip = iff(<r + 13) /2 + ip±(ao — cr 3 )/2. 

The Usadel equation is subject to appropriate bound- 
ary conditions at S/HM interfaces, which for the diffusive 
limit can be written in the form 11 
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- ( ~ l.r -l.r \ 

9,\m' ,g HM ) 



G 1 ^ 



[rh l ' r ,g] , (2) 



where g is the Green's function value at the super- 
conducting side of the appropriate S/HM interface (at 
x = Tds/2), a = +1(-1) at the left (right) S/HM in- 
terface and as stands for the conductivity of the film. 



,l,r 
J.r 



m l ' r (rpo(l + r 3 )/2 + m l ' r (T*p Q (l - r 3 )/2, where 
m l,r is the unit vector aligned with the magnetization 
direction of the left or right half metal. We assume 
that the half metals have opposite magnetization direc- 
tions, that is m r = —m l . The second term accounts 
for the different conductances of different spin direc- 
tions and Gmr ~ Gt,-\ — Gt,i- The third term ~ G^ 
gives rise to spin-dependent phase shifts of quasiparti- 
cles being reflected at the interface. Microscopically 11 
G\ r = 2(G q /S) J2JTk r ~ lWk r , where S is the junc- 
tion area and G q = e 2 jh is the quantum conductance. 
Summation over n means summation over transmission 
channels. T\f is the transmission probability for the n th 
channel and d<$n^ is the phase difference betweeen wave 
functions of spin-up and spin-down electrons, acquired 
upon reflection from S/HM interface (spin-dependent 
phase shift). Boundary conditions @ are only valid for 
small (with respect to unity) values of transparency T n 
and spin-dependent phase shift d<p n in one transmission 
channel. The value of d(f>n can be roughly estimated, 
for example, by modelling the barrier at S/HM interface 
by U a {x) = U a S(x). Then in the tunnel limit T n <C 1 
one obtains d(j) n ps vf{Ui — U^)/U^Ui. In general, the 
boundary conditions can contain another term propor- 
tional to G x accounting for spin-dependent phase shifts 
of quasiparticles upon transmission^. However, we are 



mostly interested in the tunnel limit, where this term can 
be disregarded with respect to Gq. 

9hm stands for the Green's functions at the half metal- 
lic side of the interface. Since in half metals a Fermi 
surface only exists for one of the spin orientations, the 
standard quasiclassical description is inapplicable. How- 
ever, half metals still allow for a straightforward quasi- 
classical treatment in the separate-band picture: quasi- 
particle trajectories simply exist only for one of the spin 
orientations 18 . If one chooses the quantization axis along 



the left HM magnetization, g£ M take the form 



-R,A 

9hm 



-iTTKT 3 (a + a<7 3 )/2 



(3) 



Here labels (I, r) are omitted for brevity and n = +1( — 1) 
for the retarded (advanced) Green's functions. The dis- 
tribution functions in the half metals are assumed to have 
the equilibrium form shifted by the applied voltages Vj jr . 
We suppose that V r = —V\ = V . In this case 
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(4) 



The self-consistent order parameter in the film is ex- 
pressed via the Keldysh part of the anomalous Green's 
function. We assume that the paring constant is non- 
zero only for the singlet pairing channel. Then the cor- 
responding self-consistency equation takes the form 



X f de 



^[f R (e)^(e)+f^(e)ip^)} ,(5) 



where A is the dimensionless coupling constant. 

Eqs. (p}-© constitute a full system of equations for 
solving the problem. For the case of a thin supercon- 
ducting film ds <C £s they can be solved analytically just 
analogously to the case of thin S/F bilayer 6 , mentioned 
in the introduction. Averaging retarded and advanced 
parts of Eq. (TTJ over the thickness of the film and tak- 
ing into account the boundary conditions, one can reduce 
the Usadel equation to an equation describing an uniform 
magnetic superconductor with an effective exchange en- 
ergy h e ff and a decoherence factor T. This equation can 
be easily solved. The corresponding anomalous compo- 
nents of retarded and advanced Green's function, enter- 
ing Eq. (|5|) take the form 



R,A 



7rA 



A 2 - [e + (xheff + inVy 



(6) 



where the decoherence factor T = (Gj. + 2G r MR + G l T + 
2G l MR )D/4asds physically describes the leakage of su- 
perconducting correlations from the film into the HM 
regions. This term is quite standard (including S/N sys- 
tems) and is responsible for the destroying of supercon- 
ductivity in thin films by the proximity effect. 
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The effective exchange energy h e ff = (G^, — 
G l ^)D /2asds is generated by S/HM interfaces. It is 
inversely proportional to the film width d$- As it was 
mentioned above, the boundary conditions ^ and, cor- 
respondingly, Eq. ^ are valid for d<j) n <C 1. Beyond this 
limit the effect of magnetic boundaries cannot be reduced 
to the effective exchange in the filnJ^, but there appear 
additional terms to some extent analogous to the mag- 
netic impurities. We assume that the condition d(p n <C 1 
is fulfilled. However, this does not mean that the result- 
ing exchange fields are small. In order to observe recov- 
ering of superconductivity, suppressed by the exchange 
field, one needs h e ft > A. For the film with ds J$ £s this 
condition is accomplished if (i) the left and right inter- 
faces are not identical, that is Gj, 7^ GL and (ii) the di- 
mensionless parameter h e f//A ~ G^^g/asds > 1. This 
parameter can be estimated as —{2Nf^G q / Sasds)d(j) ~ 
— (£g/lds)d(f>, where N is the number of transmission 
channels and / is the mean free path. From this esti- 
mate it is seen that for dirty superconductors G^ can 
generate large h e f / even for dej) <C 1 . 
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FIG. 1. (a) The dependence of A(T = 0) on eV. Solid line: 
h e ff = 0, dotted: h e ff = 1.25, dashed: h e ff = 2.5. All 
the quantities are measured in units of the bulk order pa- 
rameter Ab taken at l/r s / = and T = 0. Inset: popula- 
tion of majority and minority subbands for the given quasi- 
particle distribution, (b) The dependence of A on tem- 
perature. Gray solid line: Ab(T), black solid: Ao(T) = 
A(T,h eff = eV = 0), dashed: A(T,h eff = eV = 2.5), dot- 
ted: A(T,h eff = 2.5, eV = 3.1). For panels (a) and (b) 
l/r s f — 0. (c) The distribution function ip\.(e) for different 
spin relaxation rates: 1/t s / = (solid line), 0.01 (dashed) 
and 0.02 (dotted). h eff = eV = 2.5. ip t (e) = -^(-e). (d) 
The dependence of A on eV for /i e // = 2.5 and different spin 
relaxation rates (the same as in panel (c)). For all the panels 
T = 0.1. 

Now let us turn to the discussion of the distribution 
function. We neglect energy relaxation in the film, that 
is assume that the time r esc = which an electron 

spends in the film is much less than the energy relaxation 
time r e . Spin relaxation processes are also not taken into 
account. We discuss their influence below. Then it can be 
obtained from Eqs. ([U-S]) that the distribution function 



in the film takes the form 

, e + aeV , . 

ip a = tanh — — — . (7) 

It is worth noting here that the distribution function has 
such a one-step shape (in each of the spin susbbands) 
due to the fact that the leads are HM: the electrons from 
spin- up (spin-down) subband can flow only to/from the 
left (right) lead. This one-step form is very essential for 
the existence of the effect. In principle, the supercon- 
ductivity recovering can be also observed if one takes 
strong ferromagnets instead of half metals, but in this 
case the nonequilibrium distribution function inside the 
film is represented by a sum of the distribution functions 
coming from the left and right leads, weighted by fac- 
tors depending on the interface transparencies (this is a 
double-step structure). This would lead to only partial 
recovering of superconductivity, or even to the absence 
of the effect. So, in order to provide the appropriate 
distribution function in the film the resistances of the 
ferromagnets and S /F interfaces should obey quite strict 
conditions. 

Substituting Eqs. ©-© [<p is obtained making use of 
the symmetry relation tp^^(e) — — <Pi,f(— s)] into Eq. ([5]) 
we come to the following self-consistency equation 

WD f 

1 f de I 

x- J t{ r<! 



, e-eV] 

From Eq. (jSJ) it is obvious how superconductivity in 
the film is recovered under the simultaneous influence of 
the exchange field and the spin-dependent quasiparticle 
distribution. At h e ff — eV for each of the subbands we 
have practically the same situation as for the equilibrium 
non-magnetic film corresponding to eV = h e ff — 0. It is 
worth noting here that in the framework of the simplified 
weak-coupling model with constant pairing potential the 
maximal value of the Zeeman field, which does not de- 
stroy the superconductivity, equals to the cutoff energy 
loq. However, for high enough exchange fields this sim- 
plified model is inapplicable. For concreteness let us dis- 
cuss the phonon-mediated superconductivity. In this case 
our consideration fails for exchange fields of the order of 
the Debye energy and the correct calculation should be 
carried out in the framework of the particular phonon 
model. 

The resulting A as a function of applied voltage eV is 
plotted in panel (a) of Fig. [T] for different values of h e ff. 
It is clearly seen that the effect of superconductivity re- 
covering only takes place if the exchange field h e ff and 
the spin accumulation potential eV are very close to each 
other: their difference should be less than Ao Here Ao 
denotes the value of the superconducting order param- 
eter in the film at h e j f = 0. The physical reason for 
superconductivity recovering can be easily caught even 



sgn(e + h eff ) 
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without solving of a self-consistency equation, already 
on the level of consideration of the Cooper's problem of 
one electron pair. It can be shown that in the presence 
of the exchange field the lowest energy level Eq for a pair 
of electrons with pi = —p2 and opposite spins gets lower 
upon increasing the spin accumulation potential. Finally, 
Eq becomes exactly equal to its value for zero exchange 
field at eV = h e ff. It is worth to note here that in 
the absence of the exchange field the spin accumulation 
potential also destroys superconductivity at \eV\ ~ A . 
This effect has been studied in the literature as theoret- 
ically so as experimentallyi^— . 

At eV = and T, T — > self-consistency equation (j8]) 
has non-zero spatially uniform solution for h e ff < A.q. 
However, it is well-knownA that the uniform solution be- 
come metastable even earlier, at h e ff = Ao/-\/2 (Pauli 
limiting field) because of the fact that the paramagnetic 
state is more energetically favorable for higher exchange 
fields. In the considered case the paramagnetic state can- 
not be realized because the distribution function is cre- 
ated and supported by the external conditions in such a 
way that the populations of majority and minority sub- 
bands in the film remain equal. This is illustrated in 
the inset to panel (a) of Fig. [TJ The uniform supercon- 
ducting state is obviously more favorable than the nor- 
mal one due to the condensation energy. However, there 
can exist another possibility: in principle, under spin- 
dependent quasiparticle distribution an inhomogeneous 
superconducting state (analogous to LOFF-state) can oc- 
cur and be more favorable than the homogeneous one at 
some ranges of parameters. This issue is a prospect for 
future work. 

For the resonance value of eV res — h e ff the depen- 
dence of the superconducting order parameter on tem- 
perature is very similar to the original BCS one giving 
practically the same ratio 2A(T = 0)/T c , as illustated in 
panel (b) of Fig. [TJ As it is also represented in panel (b) 
of Fig. [TJ when eV deviates from eV res the temperature 
supresses the order parameter more sharply. 

Now we turn to the discussion of spin relaxation in- 



fluence on the effect. We assume spin- flip scattering 
from magnetic impurities to be the dominant spin re- 
laxation process inside the superconducting film at low 
temperatures. It can be taken into account by adding 
the corresponding self-energy term [—(l/2nT s f)aga, g] 
to the left-hand side of Eq. (TT]). Here a — [tr(l + 
Ts)/2 + <x*(l — T^)/2\pQ. As it is well-known, the in- 
fluence of the spin-flip scattering is twofold. Firstly, 
it "works" as a depairing factor destroying the coher- 
ence peaks and reducing the critical temperature of the 
superconductor—. Secondly, the spin-flip scattering in- 
fluences directly the distribution function reducing the 
difference ^ — tp^. The reduction can be roughly esti- 
mated as cpy - (p 3 / = ((ft - >fii)/(l + 2T esc /r s f). Here 
ip^ — ip 8 ^ is the difference in the presence of spin re- 
laxation processes, while <p-f — ip± is defined by Eq. ([7]) 
and r a f is the characteristic spin relaxation time. The 
results of the exact calculation of the distribution func- 
tion slightly deviate from this rough estimate, especially 
in the region of the coherence peaks. They are repre- 
sented in panel (c) of Fig. [TJ Obviously, the " damage" 
of the distribution function due to spin relaxation also 
suppresses the effect of A recovering. This suppression 
can be roughly viewed as the effective reduction of the 
coupling constant A — » A e // = A(l + r esc /r s /) i . The 
resulting influence of the spin-flip scattering on the order 
parameter is demonstrated in panel (d) of Fig. [TJ 

In summary, we have theoretically shown that cre- 
ation of spin-dependent quasiparticle distribution in a 
superconductor can fully compensate the pair-breaking 
effect of the Zccman field. Thus, superconductivity 
can be recovered for exchange fields well exceeding the 
Pauli limiting field if the spin accumulation potential 
eVj. — eVf ~ 2heff is generated in the superconductor. 
It is proposed that this effect can be experimentaly re- 
alized on the basis of voltage biased junction consisting 
of a thin superconducting film sandwiched between two 
half metals. 
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